Observability of radiation pressure shot noise in optomechanical systems 
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We present a theoretical study of an experiment designed to detect radiation pressure shot noise 
in an optomechanical system. Our model consists of a coherently driven optical cavity mode that 
is coupled to a mechanical oscillator. We examine the cross-correlation between two quadratures of 
the output field from the cavity. We determine under which circumstances radiation pressure shot 
noise can be detected by a measurement of this cross-correlation. This is done in the general case of 
nonzero detuning between the frequency of the drive and the cavity resonance frequency. We study 
the qualitative features of the different contributions to the cross-correlator and provide quantitative 
figures of merit for the relative importance of the radiation pressure shot noise contribution to other 
contributions. We also propose a modified setup of this experiment relevant to the "membrane-in- 
the-middle" geometry, which potentially can avoid the problems of static bistability and classical 
noise in the drive. 

PACS numbers: 42.50.Lc, 42.50.Wk, 42.50.Ct, 05.40.Jc 



I. INTRODUCTION 

An optomechanical system is characterized by an inter- 
action between light and mechanical motion 0, Q ■ The 
recent interest in such systems was initiated by the ef- 
forts to detect gravitational waves [3, |4[ . The field has 
now taken on a life of its own, especially since the obser- 
vation of quantum effects in mechanical systems is nearly 
within experimental reach. Besides the potential for tech- 
nological innovation, the possibility to study mechanical 
systems in the quantum regime [l], [|| could provide in- 
sight into the fundamentals of quantum mechanics [H, Q . 

The canonical optomechanical system consists of an 
optical cavity where one of the end mirrors is free to 
move 7, 8] . When light from a laser enters the cavity, the 
light exerts a force on the movable mirror. As the mirror 
moves, the cavity length changes, altering the resonance 
frequency of the cavity and thus the optical intensity in 
the cavity This in turn changes the optical force on the 
mirror, such that the optical and mechanical dynamics 
are coupled. Experimental studies of this system 0-[3| 
are now getting close to observing quantum effects due 
to the optomechanical coupling. 

Other experimental realizations of coupled optical and 
mechanical degrees of freedom have emerged over the last 
decade. One realization includes placing a delicate mem- 
brane inside an ordinary optical cavity with fixed end 
mirrors [H, 16]. This has the advantage of not having to 
combine the flexibility needed for the mechanical oscilla- 
tor with the rigidity of a high finesse cavity mirror. We 
will refer to this setup as the membrane-in-the-middlc 
geometry. Other examples include mechanical breathing 
modes in toroidal microresonators (l7l fl8| . nanobeams 
coupled to microwave resonators in superconducting cir- 
cuits Ql| d(| , optical forces on free-standing waveguides 
[2TI l22j |. and coupling to collective movement of cold 
atoms in an optical lattice [23, |24| . 

One of the major goals in the field of optomechanics is 
the observation of radiation pressure shot noise (RPSN) . 



This is the radiation pressure fluctuations experienced by 
the mechanical oscillator due to photon number fluctu- 
ations. Equivalently, RPSN is the quantum back-action 
of an optical displacement measurement [25]. One ap- 
proach to observe RPSN would be to measure a correla- 
tion between photon number fluctuations (shot noise) in 
the optical field and the position fluctuations of the me- 
chanical oscillator. The main obstacle to observing this 
correlation is that the shot noise induced fluctuations of 
the mechanical oscillator are typically very small com- 
pared to thermal fluctuations associated with the me- 
chanical damping. To detect such a correlation in the 
movable-mirror geometry, Heidmann et al. [26| proposed 
using two optical beams, one strong signal beam and one 
weak meter beam. The position fluctuations of the me- 
chanical oscillator induced by the signal beam would be 
detected with the meter beam. If the beams are exactly 
at the cavity mean resonance frequency, a correlation be- 
tween the fluctuations in the signal beam intensity and 
the phase quadrature in the meter beam will be due to 
RPSN alone. 

The idea presented by Heidmann et al. [26| has been 
tested in the classical domain [13] . R was shown exper- 
imentally that this scheme works when imposing addi- 
tional classical noise in the beam, but correlations due 
to the smaller quantum shot noise have not yet been 
observed. A more indirect observation of RPSN was 
achieved in a cold atom experiment [23] through the mea- 
surement of RPSN-induced heating of the atomic ensem- 
ble. 

In this article, we report theoretical studies of an ex- 
periment similar to the one proposed in Ref. [26|. We 
consider a two-sided cavity, where the cavity mode is 
driven by one rather than two beams and coupled to a 
mechanical oscillator. We examine the cross-correlation 
between two quadratures of the cavity output field, 
studying the general case where the beam is not at reso- 
nance with the cavity. The qualitative features of the var- 
ious contributions to the cross-correlation are explored, 
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and we determine under which circumstances the RPSN 
contribution can dominate over the contribution from 
thermal noise. This is quantified by simple figures of 
merit that can be used to determine the parameter val- 
ues needed in such an experiment. Furthermore, we 
compare the RPSN contribution to the unwanted con- 
tribution from classical noise in the laser. Finally, we 
also propose a modified setup for the membrane-in-the- 
middle geometry in which two optical modes couple to 
the mechanical oscillator with opposite signs. This has 
the potential to overcome the problem of static mechan- 
ical bistability [ill [H[ , as well as to significantly reduce 
the contribution from classical laser noise. 

The article is organized as follows: In Section HH we 
present the model and briefly discuss its properties. Sec- 
tion lllll gives details on the detection scheme, presents the 
various contributions to the cross-correlation, and shows 
how to minimize the thermal noise contribution. The 
comparison of the quantum and thermal contributions is 
given in Section llVl whereas Section IVl contains the com- 
parison to the contribution from classical laser noise. In 
Section [VTl the new two-mode setup is presented. Our 
conclusions are presented in Section [VII| and mathemat- 
ical details can be found in Appendix [A] and |B] 



II. MODEL 

We consider one optical mode in a two-sided cavity 
coupled linearly to the position of a mechanical oscillator. 
The Hamiltonian is 



scattering of photons into other optical modes or absorp- 
tion in the end mirrors or membrane. The input mode 
ain,M is the associated quantum noise operator. The total 
decay rate of the cavity mode is the sum of the contribu- 
tion from each port, i.e. k = kl + «r + km- The optical 
output modes are given by 



a ou t,i(t) = \fKi&(t) - a in j(t) , i = L,R,M 



(3) 
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FIG. 1: (color online). Schematic overview of the model. The 
optical cavity mode a is coupled to three input modes and 
to the position fluctuations z of the mechanical oscillator, 
represented in this case by a membrane in the middle. We 
imagine the cavity being coherently driven by a laser on the 
left hand side (see Figure [2]), such that the input mode ai n ,L 
is the sum of a mean amplitude, classical laser noise, and 
quantum noise. The mode fim.R simply represents quantum 
noise entering the cavity from the right hand side, whereas 
fim.M represents quantum noise associated with other types 
of decay. 



H = hoj M c i c + h(ujc + Az) (a f a - (a f a)) +H K + H 1 (I) 

where c(a) is the annihilation operator for the mechan- 
ical(optical) oscillator with frequency ojm(wc). The di- 
mensionless position operator of the mechanical oscillator 
is z = c + c? and A is a coupling constant. H K and H 1 
contain the couplings to an optical and mechanical bath, 
respectively, and describe drive and/or decay of the os- 
cillators. In the membrane- in-the- middle setup [Tol. [T^. 
when the cavity frequency depends linearly on membrane 
position, this model is valid in the limit of small mem- 
brane reflectivities. 

In the Markov approximation, input-output theory (29l . 
l3(ij gives the quantum Langevin equations 



+ iuJc ) o, — iAza + t/k~Zci 



(2) 



c = - ^~ + iuJuj c — iA (a* a — (a^a)) + y^/rj . 

Here, ai n ,L(oi nj R) is the input mode at the left(right) 
hand side of the cavity. See Figure [T] for a schematic 
picture of the input-output modes. The decay rate due 
to a finite transmission of the end mirrors is characterized 
by kl and kr. We have included a third decay channel of 
strength km, which can describe optical loss, e.g., due to 



The mechanical input mode f] describes noise from the 
mechanical bath. In the case of high mechanical quality 
factor wm/7, we can assume 



(n th + l)8(t 



•0 



(4) 



where n t h is the mean number of phonons in the absence 
of optomechanical coupling, determined by the tempera- 
ture of the mechanical bath [34| . 

We will assume that the cavity mode is driven at the 
frequency wd from the left hand side of the cavity. See 
Figure [5] for a schematic setup. Thus, we can write 



a in ,L(t) = e-^ f (%(i)+6,(t)) 



(5) 



&in,i(*) 



— iuut 



Zi(t) , i = R,M 



We write the coherent state amplitude as <Xi n (i) — ao + 
5x(t) + i5y(t), where ao is defined to be constant and 
real, such that the real functions Sx(t) and Sy(t) describe 
classical amplitude and phase fluctuations of the drive, 
respectively, e.g., arising from technical laser noise. The 
noise operators £j satisfy 



(ii(t)tl(t')) - (nc + l)*(t-f) 
(#(t)6(f)) = n c 5(t-t') 



(0) 
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to a good approximation, where nc is the thermal occu- 
pation number at the cavity resonance frequency [3(| ■ We 
will make the experimentally relevant assumption that 
hue 3> k-gT and set nc to zero. We will also assume 
that the classical amplitude and phase noise is white, 
since we are only interested in the noise around the me- 
chanical frequency ujm- We denote the strength of the 
classical noises by Cx and Cy. 



(Sx{t)Sx(t')) 
(Sy(t)Sy(t')) 



C x S(t 
C Y S(t 



t') 





(7) 



Cx and Cy are typically both functions of the laser 
power. We assume the amplitude and phase noise to 
be uncorrelated. 

In the frame rotating at the optical drive frequency, we 
write the cavity amplitude as the sum of its mean value 
and a fluctuating part, a(t) = e~ l " D *(a + d(t)). When 
we neglect the small terms <$d and dz, the linearized 
equations of motion become 



d = - ( ~ - iA ) d - iaz + V^l ( Sx + iSy + £l ) (8) 



+ 



/ kr£,r + \/km£m 



Here, the detuning A = wd — wc, the coupling a = Aa, 
and the mean cavity amplitude a = (a) = ^fn^a® / (k / 2 — 
iA). Note that at nonzero detuning A, the mean cav- 
ity amplitude a is phase shifted relative to the incident 
amplitude by a phase (f> given by 



= arctan ■ 



2A 



(9) 



The solutions to the equations of motion are given 
in Equation (|A1|) of Appendix [A] They can be ex- 
pressed in terms of the mechanical susceptibility [35j 
Xm[w] = [7/2 — i (w - um)] - i the cavity susceptibil- 
ity XcM = [n/2 — i{uj + A)] -1 , and the optomechanical 
"self-energy" 



= -i\a\ 2 (xcH - XchH) 



(10) 



We will let 7 opt denote the additional damping of the 
mechanical oscillator due to the optomechanical interac- 
tion. In the weak-coupling limit where |7opt| <S k,wm, 
the mechanical frequency and damping rate are shifted 
by (Jcjm = ReS[wM] and 7 op t = -2IhiE[wm], respec- 
tively [31|. We will assume to be in this limit and let 



u M — ui M + Sujm 
7 = 7 + 7o P t 



(11) 



denote the effective mechanical frequency and damping 
rate. The effective mean phonon number becomes 



7 nth + 7o P t ^o P t 
7 



(12) 



where n op t is given in Equation (|A14|) . In the absence 

of classical laser noise, n opt — — (4wmA|xc[wm] | 2 ) 
is a measure of the effective temperature of the RPSN 
[3fl [3l1 |. At a sufficiently large positive A, the effective 
damping 7 becomes ne gati ve and the mechanical oscil- 
lator becomes unstable 36]. At negative A, the optical 
damping 7 opt is positive, leading to cooling of the me- 
chanical motion. 



III. DETECTING RADIATION PRESSURE 
SHOT NOISE 

As proposed by Heidmann et al. [26J, RPSN on the me- 
chanical oscillator can be detected by a cross-correlation 
measurement of the outgoing beams from the cavity. We 
consider a similar detection scheme here. The general 
idea is presented in Section fill A| and Section QlIB] pro- 
vides a motivation. In Section IIII CI we show why me- 
chanical thermal noise poses a problem and how it can be 
avoided. In Section IIIIDl we present a modified version 
of the experiment, and give general expressions for the 
various contributions to the cross-correlation function. 



A. General scheme 

Figure [2] shows a simplified schematic of the most gen- 
eral setup we consider. 
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FIG. 2: (color online). Schematic overview of the experiment 
considered. The cavity mode is coherently driven by a laser 
from the left hand side. One quadrature Xe R of the transmit- 
ted beam and one quadrature Yg L of the reflected beam are 
measured through homodyne detection. The RPSN felt by 
the mechanical oscillator can be detected through the cross- 
correlation of the fluctuations in the two quadratures. 



We assume that the quadrature 

X 9R (t) =e i ^ t -fe)a QUt , R (t)+e-<^*-^)aL t , R (t) (13) 

of the transmitted beam is measured through homodyne 
detection, with some as yet unspecified local oscillator 
phase The fluctuations around the mean value of 
X@ R (t) can be written 



SXg R (t) 



<W(i) + e l9R dl ut At) , (14) 
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where d ou t,R,(t) = v / KRii(i) — £r(£). Similarly, on the left 
hand side of the cavity, we assume that the quadrature 
fluctuation 



-iO L 



^out,L 



(t)+e 



'dint At) > 



(15) 
Here, 



is measured through homodyne detection. 
4ut, L (t) = ^d(t) - (5x(t) + i5y{t) + &(*))■ 

For the detection of RPSN, we examine the Fourier 
transform of the cross-correlation of the two quadrature 
fluctuations 8Xg R and 8Yg^ defined above. Since they in 
general do not commute, it is important to note that a 
measurement will correspond to the symmetrized cross- 
correlation. Hence, we study the function 

i r°° 

SH = g / dte^{{8X 6H {t),SYeM}) (16) 



1 

47T 



dJ {{8X 9R [u>],SY eij [u)']}) 



where the brackets denote the anticommutator. In prac- 
tice, S[lo] can easily be measured by multiplying the com- 
plex Fourier-transforms of time traces of 8Xg R (t) and 
SYg^t) d?) ■ We choose not to normalize S[w] by the 
auto-correlators as in Ref. [26], since we are not only 
looking for a nonzero value, but are also interested in the 
qualitative frequency dependence, which is complicated 
by normalization. Also, comparing the size of the signal 
for different parameter values can be of interest when it 
comes to the problem of technical noise. 



B. Motivation 

The interaction part of the Hamiltonian ([I} can in the 
linearized case be written H{ nt — —IqzF, where the opti- 
cal force operator is 



F = — 



(17) 



and Iq is the size of the zero point fluctuations of the 
oscillator. We now briefly motivate how the detection 
of optomechanical correlations, i.e., correlations between 
the optical force F and the mechanical oscillator position 
I, can be inferred from the measurement of S[ui]. 

We consider a simplified situation where the laser is 
on resonance with the cavity (A = 0). In addition, we 
neglect the quantum vacuum noise £i and the classical 
phase noise 8y{t). Let us at first assume that the clas- 
sical amplitude noise 5x is simply a periodic signal at 
frequency wn, i.e., 8x(t) ~ cos(wn0- Furthermore, we 
let #r = 0, such that Xg R is the intensity quadrature on 
the right hand side, and #l = tt/2, such that Fg L is the 
phase quadrature on the left hand side. Finally, we as- 
sume that the cavity decay rate is much larger than the 
oscillation frequency of 8x and the resonance frequency 
of the mechanical oscillator, i.e., k 3> lo^,ujm- We drop 



hats on operators in the discussion below, since all quan- 
tum effects are neglected. 

In this simplified case, the quadrature 8Xo is simply 
proportional to the optical force F(t), which again is pro- 
portional to 8x(t). On the other hand, the reflected phase 
quadrature is proportional to the position of the mechan- 
ical oscillator, 8Y^/ 2 ~ z (t)- The cross-correlation be- 
tween the two quadratures becomes 



({6X (t),5Y„ /2 (Q)}) ~ (F(t)z(fi)) 



(18) 



This means that the function S[ui] is the Fourier trans- 
form of the cross-correlation between the optical force 
and the position of the oscillator. 

We can take this simple analysis one step further and 
consider what kind of qualitative behaviour we can ex- 
pect for the function S[u>}. The force is simply propor- 
tional to the signal 8x, such that F(t) ~ cos(wNi)- If we 
neglect the mechanical bath and assume 7, |wn — wm| ^ 
ujm, the position of the mechanical oscillator is 



~ Ixm[wn]| cos (w N t - A(w N )) 



(19) 



This is phase shifted relative to the force, where A(wn) 
is given by 



A(wn) = arctan ■ 



7/2 



ujm — Wn 



(20) 



As expected for a damped and driven harmonic oscillator, 
A w for frequencies ^> 7, A — n/2 for cj^ = ^m, 

and A w ir for ojn — 3> 7. 

Interpreting the expectation value in Equation (fT6)) as 
a time average, one arrives at S[w] — S[u] + S*[— u], with 



S[u 



ojm -wn + 17/2 
(7/2) 2 + (w N - com) 2 



8(lu — cjn) 



(21) 



Remember that we assumed 8x(t) ~ cos(uj-^t). In the 
relevant situation where 8x(t) represents white noise, 
the force is driving the oscillator at all frequencies 
simultaneously. The real part of the Fourier trans- 
formed cross-correlation S[uj] will then be proportional to 
|XmMI cosA(w), which has a sign change at the mechan- 
ical frequency wm- The imaginary part will be propor- 
tional to |xm[w]| sin A(w), which is a Lorentzian centered 
at uj M - 

The above analysis gives an idea of how optomechani- 
cal correlations manifest themselves in the function S[ui]. 
However, we only considered a very simplified case, and 
the general situation is more complicated. In addition, 
we only took classical noise into account, but for the pur- 
pose of detecting RPSN, we are looking for correlations 
between quantum optical noise and oscillator position. It 
turns out that in the special case we analyzed above, the 
result is valid also for quantum noise. However, we will 
see that optomechanical correlations of quantum origin 
can actually be distinguished from their classical coun- 
terparts under certain conditions. 
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In general, to affirm that photon shot noise in the cav- 
ity is influencing the mechanical oscillator fluctuations, it 
is simply enough to prove that correlations exist between 
the position operator z and the vacuum noise operators ^ 
which are the sources of the shot noise. We will see that, 
in the general case, S[uj\ contains such correlation func- 
tions, and we will determine under which circumstances 
additional terms can be neglected. 



C. The problem of thermal noise 

In the simplified example in the previous section, the 
transmitted intensity quadrature was independent of the 
oscillator position. This is always the case when the de- 
tuning A is exactly zero, as noted by Hcidmann et al. 
(26[. In general however, both quadratures SXg R and 
SYg L depend on the oscillator position z. This gives rise 
to a term in S[uj], denoted S z ,zM below, that is propor- 
tional to the spectral density of the mechanical oscilla- 
tor, which is typically dominated by thermal noise due 
to the mechanical bath. If this term is much larger than 
the terms originating from optomechanical correlations, 
RPSN detection by this method becomes difficult. How- 
ever, we now show that the thermal noise contribution 
can be made to vanish not only for zero detuning, but 
for any detuning A. This will permit the observation of 
the optomechanical correlations. 

From Equation ([5]) , we see that the cavity field fluctu- 
ations d depend on the position operator z, the quantum 
vacuum noise £j, and the classical laser noise Sx, Sy. Since 
the equations of motion are linear, we can focus on the 
part that depends on z and ignore the other terms. The 
cavity field fluctuations then become 



d(t) 



d re -(»/2-iA)(i-T) |( - T ) ( 22 ) 



Let us first consider the bad cavity limit n 3> wm, where 
the cavity field follows the motion of the oscillator adia- 
batically. In this limit, a general quadrature of the cavity 
field becomes 



-if) 



d{i) + e l9 d\t) 



2|a|sin(0-2(/>) 



m . (23) 



This shows that in the bad cavity limit, choosing 0r(0l) 
to be 2<j) makes Xg R (5Yg h ) independent of the position z, 
such that the thermal noise contribution in S[u>] vanishes. 
This also connects smoothly with the resonant case A = 
0, where the intensity quadrature (0 = (f> = 0) becomes 
independent of z. 

If the experiment is not in the bad cavity limit, the 
picture is more complicated, but it is still possible to 
avoid the thermal noise. We show this by first assuming 
that the motion of the oscillator is given simply by z(t) ~ 
cos(a>Nt), i.e. , a harmonic oscillation at frequency wn- 
Again, we drop operator hats as everything is considered 



to be classical. The cavity quadrature given by becomes 

e- ie d(t) + e i6 d\t) = -2\a\\f(e)\z{t-p(9)/us) (24) 

where f(9) = Xc[^n]Xc i~ w n]((^/2 — iu>w) sin0 — A cos 9) 
and 9 — 9 ~ <fr is the quadrature phase measured relative 
to the phase of the intensity quadrature. We observe that 
there is a phase shift p(9) between the oscillator motion 
and the cavity field fluctuations, defined by cxp(ip(0)) — 
f(9)/\f(9)\. Note again that in the bad cavity limit k ^> 
wn, the choice 9 = 2<j> makes the prefactor (1/(0)1) vanish. 

In the general case, when examining the cross- 
correlation (8Xe R (t)6Yg^(0)), one finds that the term 
symmetric in time, which corresponds to the real part 
of S[ui], vanishes when p(0r) — p(6h) = ±7r/2, whereas 
the antisymmetric part, corresponding to the imaginary 
part of S[uj], vanishes when p(0r) — p(&l) = 0, tt. The 
latter is always the case when 0r = 0l- The former cri- 
terion demands that the two quadrature phases satisfy 



A 



cos(0r - 9 L ) ~ kA sin(0R + L ) 
cos(0r + L ) = . (25) 



By choosing one of the quadratures, this equation gives 
the other quadrature for which Re5[aj] will vanish. 
The physical interpretation is that the two quadratures 
8Xq r (t) and SYg R (t) then measure orthogonal quadra- 
tures of the mechanical oscillator's motion. 

In the toy example above, we made the assumption 
z(t) ~ cos(wn^)- I n reality, the oscillator motion is a 
noisy signal and not restricted to one frequency. If we 
still want the real part of the cross-correlation S[cj] to 
vanish, Equation (1231) must be fulfilled for all frequencies 
wn. This is only possible if either 0r or 0l is equal to 0, 
i.e., the intensity quadrature in the cavity. The natural 
choice is then 0r = <j>, since this can be achieved simply 
by replacing the homodyne detection on the right hand 
side with a photomultiplier and recording the intensity 
fluctuations. For the remainder of Equation ([^5)1 

then dictates the other quadrature to be 0l = 2(f>. To 
locate the correct 0l, one possibility is to drive the os- 
cillator mechanically and look for the quadrature phase 
where the strong signal from the driven oscillator disap- 
pears in Re S[lj]. This null can be maintained via a servo 
loop throughout the course of the experiment provided 
that the external mechanical drive is at a frequency well 
away from the frequencies of interest. 

Let us also mention that for a mechanical oscillator 
with a high quality factor wm/t, the requirement that 
one of the quadrature phases must be <f> can possibly 
be relaxed by only demanding Equation (|25p to be valid 
for cjn = <^m- In that case, the thermal contribution 
to Res' [a;] only approximately vanishes for frequencies 
close to the mechanical frequency. The calibration pro- 
cedure mentioned above should also work in that case, 
but then only for drive frequencies close to the mechan- 
ical resonance frequency. Using 0r = <f> and 0l = 2cf> 
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does however also have the advantage that both angles 
are known, which is useful when trying to fit theoretical 
expressions to experimental results. 



D. General expressions 

From the above discussion of the thermal noise contri- 
bution, we can conclude that it is favorable to let SXg R be 
the fluctuations in the intensity quadrature, both for zero 
and nonzero detuning A. This means that Or = <p, and 
we will restrict ourselves to this situation from now on. 
We also rename the quadratures by letting 8X$ — > SX 
and 8Yg^ — > 8Yg. Hence, we consider a modified version 
of the experiment, shown in Figure [3l where the inten- 
sity fluctuations 8X are measured by a photodetector. 
Let us however mention that the general expressions for 
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FIG. 3: (color online). Modified version of the experiment 
considered, where the transmitted intensity fluctuations are 
measured by a photodetector. 



the contributions to S[uo] presented below are valid also 
for arbitrary Or, but with different coefficients. 

In terms of z, £j, 8x and 8y, the correlation function 
(8X[uj]8Yg[uj']) can be written as a sum of five contribu- 
tions: 



(SJt[u]SY e [u/]) = (5XMJy e [u/]) q , q 

+ (SX[u]SY [u J '}) chcl + (SX[u]SY [u/]) v 
+ (SX[u]SY [u J '}) chz + {8X[u]8Y e [J])^ . 



(26) 



Here, the label "q" refers to the quantum fields £j, the 
label "cl" refers to the classical fields 8x, 8y, and "z" 
refers to the position i of the mechanical oscillator. We 
have grouped the terms in this way to identify the sig- 
nature of RPSN. The first term contains correlators of 
the kind [a/]), with i = L,R,M. It turns out 



that (SX[u]5Y [uj'}) c 



0. Thus, in the absence of 



the nonlincarity introduced by the radiation-pressure- 
induced motion of the mechanical oscillator, there is 
no correlation between quantum noise in the outgoing 
quadratures on each side. The second term is the cor- 
relation between the classical laser noise in the two 
quadratures, containing the correlators (Jx[w]<5x[a;']) and 
(6y[uj}5y[uj'}} . The third term, (8X[u}]8Yg[u}'}}^ z , con- 
tains correlations between the quantum vacuum noise 



and the oscillator position, (^i[uj]z[uj']) , which are the op- 
tomechanical correlations that we would like to detect. 
The term (8X[ui]SYg[uj'}) c i. z contains correlation func- 
tions of the kind ((5x[w]z[w']), which are also optomechan- 
ical correlations, but due to classical noise in the drive, 
not photon shot noise. The last term (~ (z[u;]z[u/])), 
discussed in the previous section, is proportional to the 
position spectral density and appears as a result of de- 
tecting the mechanical oscillator position fluctuations in 
both quadratures. Although it vanishes at zero detuning 
A, this term is generically non-zero. 

By using Equations ([J}, ©, © and (|A1|) . we can 
evaluate the expectation values and write the cross- 
correlation S[uj] as 



S[oj] = Sq.zM + ScLclM + S c l tZ [u)] + s ZtZ [u] 



(27) 



where the terms correspond to the terms in Equa- 
tion (|26[) . The correlation function S[u>] is a complex 
quantity, where the real(imaginary) part is symmet- 
ric (antisymmetric) in lj. We will later use the abbre- 
viations 



R[u] 
I[ U ] 



ReS[ui} 
Im S[oj] 



(28) 



and similarly for the various contributions in Equation 
(|271) . i.e., i?q,zM = ReS'q^H etc. 

The general expression for S[ui] can be found in Ap- 
pendix |XJ We are only interested in its behaviour around 
the mechanical frequency com , where the oscillator is most 
susceptible to the radiation pressure noise. Considering 
a mechanical oscillator with a high quality factor, we 
assume 7 <C k, ujm and restrict ourselves to frequencies 
where \u> — ljm\ *C This greatly simplifies the ex- 

pressions. The (z,z)-contribution becomes 



S'z. 



(th) 



il 



(th) 



{uo-Cjm? + (7/2) 2 



(29) 



The real constants and l[ th ^ are given in Equation 

(|A11[) of Appendix [AJ This means that in the vicinity of 
the mechanical frequency ujm, both the real and imagi- 
nary parts of S^zM are Lorentzians centered around cDm 
with width 7. This is as expected, since S^zM origi- 
nates from the spectral density of the mechanical oscil- 
lator. This term contains the thermal noise contribution 
to the cross-correlator S[oj], as the spectral density of the 
oscillator is typically dominated by thermal fluctuations. 
The term we would like to detect is 



Sq,zM 



4 q) + il™ + R™ [lj - u M ) 
(w - Cj m ) 2 + ( 7 /2) 2 



r(q) 



(q) 



where the real constants R^' , r!% , can be found in 
Equation (|A15|) . The imaginary part of Sq jZ [u;] is also 
a Lorentzian. The real part is a sum of two terms, one 
proportional to the above-mentioned Lorentzian and one 



(30) 
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term which changes sign at Qm ■ We note at this point the 
important fact that when [R^/R^ | is small compared 
to 7, the real parts of Sq, z [w] and 5 Z , Z M are qualitatively 
different and in principle distinguishable. The term con- 
taining correlations between classical noise in the drive 
and oscillator position becomes 

It is similar to S qtZ [uj], but with different coefhcients and 
an additional term in the imaginary part. An important 
reason for the difference between the quantum and classi- 
cal contributions iS qjZ [u;] and Sd.zM is that classical noise 
enters the cavity only from the left hand side, whereas 
quantum noise enters through all ports. Thus, an asym- 
metric cavity with kr > kl can be favorable in terms of 
increasing the relative importance of the quantum ver- 
sus classical contributions. Apart from this, there is no 
strong dependence on the relative size of Klj «r and km- 
All plots in this article therefore refer to the case % = 
and «l = kr = k/2. The term 5 c i, c iM has no sharp 
features around the mechanical frequency and produces 
only a smooth background. 

The detection of RPSN now comes down to being able 
to identify the presence of the term S'qzjw] in the total 
signal S[us] being measured. 

IV. QUANTUM VERSUS THERMAL 
CONTRIBUTION 

In this section, we ignore classical noise in the drive, 
such that 5 c i,zM = S^clM — 0> and focus on the two 
remaining contributions. We have observed that i? z z [w] 
and i?q jZ [w] in principle can be distuingished even at 
A ^ 0, whereas the imaginary parts I ZtZ [ui] and / q , z [w] 
are indistinguishable. However, i? Z;Z [w] is proportional to 
the mean phonon number tsm, which, depending on the 
temperature T of the mechanical bath, can be a macro- 
scopic number. The contribution i? qiZ [w] will therefore 
be negligible in most cases. In other words, the fluctu- 
ations of the oscillator due to the mechanical bath are 
usually much larger than those caused by a fluctuating 
photon number in the cavity, and consequently the latter 
is negligible. However, we are not measuring the position 
fluctuations directly, but rather the cross-correlations be- 
tween the outgoing quadratures SX and SYg. In Section 
IIII Ci we observed that at special points in parameter 
space the term i? z , z [w] vanishes. Heidmann et al. (26j 
noted that this occurs at resonance, i.e., A = 0, where 
8X is independent of z. We note that for non-zero de- 
tuning, it also occurs at the critical angle 9 = 9 Cl where 
9 C = 2<p modulo 7r and given by 

A/k 

9 C — arctan ^ + kn , k g Z , (32) 

1/4 -(A/k) 2 

This equation can be viewed in two ways. It gives the 
appropriate detuning A to make i? ZjZ [w] vanish for a given 



cavity linewidth k and angle 9. For example, when 9 = 
7r/2, this detuning is |A| = k/2. Equivalently, Equation 
(f3"2")l tells us which quadrature to measure on the left 
hand side of the cavity for a given A/k. Figure [?] shows 
the critical angle 9 C as a function of detuning A. 




A/k 



FIG. 4: (color online). The critical angle 6 C as a function of 
detuning A. The angle 9 C is multivalued, since 6 — !> 6 + n 
gives 5Ye(t) — > —5Y g (t) and hence S[u>] —> —S[u]. 

Although the criteria for vanishing i? ZjZ [w] can be de- 
termined quite easily, it does not mean that experimental 
detection of R qyZ [ui] becomes straightforward. There is al- 
ways an uncertainty in experimental parameters as well 
as drift and fluctuations in laser frequency. The question 
is whether one can get sufficiently close to these criteria 
in order to claim that the correlation R^ z [ui] has been de- 
tected. In the following, we therefore aim to determine 
what sufficiently close means in quantitative terms. 

A. Zero detuning 

The unwanted contribution i? ZjZ [w] changes sign with 
A and vanishes when A = 0. One might therefore be able 
to make R ZtZ [oj] small enough that i?q !Z [w] is the domi- 
nant contribution to the cross-correlation measurement 
by choosing the detuning as close to A = as possible. 
Of course, a fluctuating and/or drifting A might pose a 
challenge in an actual experiment. 

In the limit |A| < k,cj m , the ratio |i?i q) /(#2 q) 7)l -> 0, 
such that the quantum contribution i? qjZ [u;] will have a 
sign change at (Dm- In fact, 5 q . z takes exactly the form of 
the simple force-position correlation discussed in Section 
IIII B[ but with um — > wm and 7—^7. Also, 5 qiZ [o;] is 
proportional to sin# in this limit, such that 9 = 7r/2 will 
maximize the signal, as was also noted in Ref. [26|. 

Figure [S] shows the two contributions i? qjZ [o;] and 
i? z>z [w], as well as the total signal R[oj] = R q _ z [uj] + 
i? ZiZ [w], for a detuning A = —0.01k. The parameters cho- 
sen are relevant to the present membrane-in-the-middlc 
setup [UGl. 

In the right panel, the mean number of 
photons riphoton = |a| 2 in the cavity is 100 times larger 



8 



than in the left panel, corresponding to a difference in 
laser power. The abscissa shows the frequency devia- 
tion u> — ojm in units of the bare mechanical damping 7. 
In the case displayed in Figure El we are in the regime 
where 7 op t,^M 7, where Slum is the "optical spring" 
frequency shift and the linewidth increase 7 op t is a result 
of optical cooling of the mechanical motion associated 
with the detuning being negative [31j . In the top panel, 
i?q jZ [w] is shown, and we observe that it is dominated 

by the R% -term in Equation (|30|) . In the middle panel, 
the Lorentzian i? ZjZ [w] is shown. Increasing the photon 
number by a factor of 100, i.e., going from the left to 



the right panel, we see that both R, 



q,zL l 



and R 7 , 



be- 



come wider and are shifted further from ujm due to the 
increase in 7 opt and Sujm- We observe that the peak-to- 
peak value of R q ^[uj] is largely unchanged, whereas the 
height of R ZtZ [uj] is smaller when the photon number is 
increased. The latter is due to increased optical cooling 
of the mechanical oscillator. In the lower panel, the total 
signal R[cj] is shown. On the left hand side, the con- 
tribution i? ZjZ [w] is dominant and the total signal looks 
almost like a Lorentzian, albeit with some asymmetry. 
On the right hand side, R qyZ [ui] dominates, and the total 
signal has a sign change and is clearly not a Lorentzian. 
Assuming there is no classical noise in the drive, a signal 
R[w] as in the lower right panel of Figure [5] is a signature 
of correlations between photon shot noise in the cavity 
and position fluctuations of the mechanical oscillator. 

From Equations (f5T7|) and (l2l))) . we can derive the peak- 
to-peak value P q of the contribution i? q z [u;] and the 
height M of the Lorentzian i? z . z [u;]. To be able to detect 
correlations between photon shot noise and oscillator po- 
sition, an important figure of merit is the ratio P q /AI. 
We find that for |A| <§C k,uju, 



1 q 
M 



(n/2) 2 



2(n M + 1/2) «]A| 



(33) 



where the effective phonon number um is given in Equa- 
tion (|T2"j). When this ratio becomes larger than one, the 
total signal R[oj] is dominated by R, 



In Figure [5j 

this ratio is 0.34 (left panel) and 5.42 (right panel). 

Let us assume that the parameters k>M, K, 7, nth are 
fixed and we are free to change the detuning A and the 
optomechanical coupling |a|. The latter can be varied 
by changing the input power and hence the intracavity 
photon number. First, consider the regime 7 op t <C 7, 
where nyi ~ n t h- In that case P q /M is simply inversely 
proportional to |A|. In the opposite regime, 7 op t 3> 7, 
Equation (|3"3")l becomes 



1 q 
M 



2w M |o:| : 



7™th 



(k/2) J 



(34) 



when assuming um ~> 1- Note that in this regime (|A| <C 
k,um, 7 *C 7o P t) 7 the figure of merit is independent of 
A. The reason is that an increase of A is compensated 
by additional cooling of the oscillator. The ratio can be 
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FIG. 5: (color online). The contributions i? q . z [u] (upper 
panel) and i? ZlZ [ij] (middle panel), and their sum R[lu] (lower 
panel). We have used A = — 0.01k, wm/k = 3.2, uim/i = 10 6 , 
A/um = 1-15 ■ 10~ 6 and 6 = tt/2. The temperature is 4.2 K. 
The mean number of photons in the cavity (n p hoton = ] o.| 2 ) 
is 10 10 (left panel) and 10 12 (right panel). In the first case, 
the total signal is dominated by R z , z [u)], whereas in the second 
case R q , z [u}] dominates. In the absence of classical laser noise, 
a measured signal as in the lower right panel is a signature of 
RPSN. 



made larger by increasing \a\ [3g|, i.e., input power, and 
has a maximal value of 2wm/k. It is worth noting that 
a decrease of |A| in this regime does not necessarily pay 
off in an increased P q /M. For example, a decrease of |A| 
by a factor of 10 does not significantly change this ratio 
in the examples shown in Figure [5] However, the overall 
signal gets larger and one might become less sensitive to 
technical noise. In addition, when |A| gets small enough, 
we eventually reach the point where 7 op t and 7 become 
comparable and Equation (I34]) loses validity. 



B. Finite detuning 

At finite detuning A, the signal i? Z)Z [w] vanishes when 
9 = 9 C . We now investigate whether i? q , z [u;] is sufficiently 
distinguishable from i? Z)Z [w] also in this case, and how 
small |6> — 9 C \ must be for i? qjZ [u;] to dominate. 

We noted above that the correlation we seek to 
detect, i? q , z [u;], is distinguishable from i? ZjZ [u;] when 

l^i q) /(^2 q) 7)l < 1- For ~ 6c, this requirement be- 
comes 



< 1 



(35) 



M 
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This can be satisfied in several ways, for example in the 
resolved sideband limit % 3> re with |A| <C In an 
experiment where the oscillator motion is dominated by 
thermal fluctuations, the quadrature 9 = 9 C can be lo- 
cated by looking for an overall sign change in the total 
signal R[u>], either by varying 8 or, if 9 is fixed, by varying 
A. Exactly at 9 = 9 C , one should be able to see the asym- 
metric features originating from i? qjZ [w]. As mentioned 
earlier, the calibration process for locating 9 C can be im- 
proved by strongly driving the oscillator mechanically to 
temporarily amplify the (z,z)-contribution. 
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FIG. 6: (color online). The contributions P q , z [tj] (upper 
panel) and P z , z [u] (middle panel), and their sum R[to] (lower 
panel). We have used A = — 0.5k, ljm/k = 3.2, ujm /j = 10 6 
and A/lom — 1.15- 10 -6 . The temperature is 4.2 K. The mean 
number of photons in the cavity (n p i loton = |a| 2 ) is 10 10 . The 
angle 6 is 6 C — 7r/10 (left panel) and 8 C — 7r/100 (right panel). 
In the first case, the total signal is dominated by R z , z [lu], 
whereas in the second case, P q , z [u;] dominates. 

In Figure [51 we show the real part of the cross- 
correlation R{uj] and its contributions for A = — re/2. In 
this case, 9 C = ir/2. In the left panel 9 = 9 C — 7r/10, and 
in the right panel 9 = 9 C — 7r/100. We see that the un- 
wanted contribution i? z . z [u;] dominates in the first case, 
whereas R q _ z [uj] is dominant when 9 gets closer to 9 C . 

Again, the important figure of merit is the ratio be- 
tween the peak-to-peak value P q of the contribution 
-Rq iZ [w] and the height M of the Lorentzian i? ZjZ [w]. For 
56 = \9-e c \ -C 1, it is 

M 4(n M + l/2)[(re/2) 2 + A 2 ]^ ' [ ' 

In Figure [51 this ratio is 0.54 (left panel) and 5.44 (right 
panel) [39j. 



Figure [7] shows the ratio (l36l) as a function of A at a 
temperature of 300 K with 59 = 7r/100. We have used 
um/k = 1.6 and wm/k = 6.4 in the left and right panel, 
respectively. Figure [3] shows the same ratio, only at T = 
4.2 K, giving significantly higher values. We observe that 
for small um/k, the ratio is not very dependent on A, 
but for larger values of lum/k it has a peak at A = —re/2. 
For high temperatures, it also has a peak at A = -wm. 
This is in accordance with the fact that cooling of the 
mechanical motion is most effective at A = -wm when 
wm/k is large [U H2[. Note however that the condition 
A = -wm with wm S> k does not satisfy Equation (|35[) . 
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FIG. 7: (color online). The ratio between the peak-to-peak 
value P q of i? qiZ [aj] and the height M of 7? z , z [w] for ojm/k ~ 1.6 
(left panel) and lum/k ~ 6.4 (right panel). The deviation from 
the critical angle 9 C is 89 — 7r/100, um/j = 10 6 , A/ujm = 
1.15 • 10 -6 , n p hoton = 10 10 , and the temperature is 300 K. 
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FIG. 8: (color online). The ratio between the peak-to-peak 
value P q of Pq iZ [oj] and the height M of 7? z , z [w] for cjm/k = 1.6 
(left panel) and um/k = 6.4 (right panel). The deviation from 
the critical angle 9 C is 89 — 7r/100, um/j = 10 6 , A/ujm = 
1.15 ■ 10 -6 , n p hoton = 10 10 , and the temperature is 4.2 K. 



From the above, we can conclude that it might be pos- 
sible to observe the correlation between photon shot noise 
and oscillator position using a negative detuning A and 
9 = 9 C , especially in the resolved sideband regime where 
Wm/k is large. We should point out that 9 C is suscep- 
tible to fluctuations in the detuning A, as can be seen 
from Equation (|32p . From this point of view, it would 
be beneficial to choose |A| large compared to the cav- 
ity linewidth re, since d9 c /dA —> when |A|/re becomes 
large (see Figure BJ. However, the choice of the detuning 
A is also limited by other requirements, such as Equation 

The potential insensitivity to laser frequency fluctu- 
ations for large |A|/re could offer an advantage of this 
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method versus performing the experiment at resonance. 
Another advantage is a lower effective phonon number 
riM due to greater optical cooling of the mechanical mo- 
tion. On the other hand, the reduced overall size of the 
signal, which reduces the signal to (technical) noise ratio, 
might be a disadvantage. 



QUANTUM VERSUS CLASSICAL 
CONTRIBUTION 



Even if a signal as in the lower right panel of Fig- 
ures [5] and [6] is detected, one needs to make sure that 
it represents a correlation between the mechanical oscil- 
lator position and photon shot noise, not classical laser 
noise. We now ignore the contribution 5 ZjZ [w] and look 
at whether or not we can distinguish the classical and 
quantum contributions. We again focus on the two cases 
where the thermal signal i? ZjZ [w] vanishes, A ~ and, 
for A ^ 0, 6 = 6 C . 




A. Zero detuning 

We saw earlier that i? qiZ [w] has a sign change at com 

when | A| < k, uj m , since \R[ q) /(R 2 q) j)\ in that limit. 
The corresponding ratio for the classical signal i? c i. z [w] 

becomes |-Ri Cl V (^2° 1 ' > t)I — w m/k when |A| <C k,wm. 
This means that in the good cavity limit wm/k S> 1 and 
close to resonance, i? c i,z[w] is approximately a Lorentzian 
(see Equation l|3"Tj) h and the classical and quantum con- 
tributions are distinguishable. This is a consequence 
of the cavity being two-sided. In the bad cavity limit 
wm/k <C 1, the classical and quantum signals i? qjZ [w] 
and i? c i [Z [w] cannot be distinguished, although it might 
still be possible to exclude the classical laser noise con- 
tribution from other types of measurement. 

The ratio between the peak-to-peak values of i? qjZ [w] 
and i?, c i jZ [w] is 

P q y/^Y+wj, 

Pel " 4k l Cx 

for |A| <C k,ojm- As motivated above, we observe that 
an asymmetric cavity with % > kl will reduce the rela- 
tive contribution from classical noise in this case. Equa- 
tion (|37l) also suggests operating in the resolved sideband 
regime, w m /k » 1. 

The two points noted above are illustrated in Figure 
[HI which shows the quantum contribution i? q z [o;], the 
classical contribution i? c i[w] = i? c i iZ [w] + i? c i jC i[a;], and 
their sum. We have used wm/k = 0.64 and ujm/k = 6.4 
in the left and right panels, respectively. The strength 
of the classical noise was chosen to be Cx = Cy = 1, 
making it comparable to the quantum noise. In the left 
panel, we observe that the classical and quantum signals 
are of the same order, and it is hard to distinguish them. 
In the right panel, the classical contribution looks more 



FIG. 9: (color online). The contributions i? qiZ [u] (upper 
panel), i? c i[u;] (middle panel), and their sum (lower panel). 
The detuning is A = —0.01k, cum/'y — 10 6 , A /ujm = 
1.15 • 10~ 6 , 6 = 7r/2 and n phot on = 10 10 . The ratio between 
mechanical frequency and cavity linewidth is wm/k = 0.64 
(left panel) and ojm/k = 6.4 (right panel). Note the quali- 
tative difference between the two contributions in the latter 
case. 



like a (negative) Lorentzian, whose height is smaller than 
the peak-to-peak value of the quantum contribution. 

We should also mention that the imaginary parts 
lq )Z [w] and .Tci,zM are qualitatively different, since the 
former is simply a Lorentzian. Examining the imagi- 
nary part I[lo] could therefore provide a way of deciding 
whether the observed signal R[uj] is of quantum or clas- 
sical nature. See Appendix I A II for details. 



B. Finite detuning 

For angles 9 close to 9 C , we know that the quantum sig- 
nal i? q , z [w] is distinguishable from R Z:Z [uj] when Equation 

([55| is satisfied. In this case, the coefficients and 

R 2 are so complicated that an expression like Equation 
([57)1 is not very helpful. However, the relative size of 
the quantum and classical contributions can always be 
examined by plotting the two contributions. Figure [TU] 
shows two examples, where A = —k/2, 9 — C = ir/2, 
Cx = Cy = 1, and where wm/k = 3.2 (left panel) and 
Wm/k = 6.4 (right panel). 

We observe that the quantum contribution dominates 
in both cases, and that a large ratio um/k seems to be 
favorable also here. 
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FIG. 10: (color online). The contributions R qtZ [u] (upper 
panel), Rc\[lj] (middle panel), and their sum (lower panel). 
The detuning is A = -«/2, lum/j = 10 6 , A/um = 1.15- 10 -8 , 
^photon = 10 10 and 6 = 6 C . The ratio between mechanical 
frequency and cavity linewidth is um/k = 3.2 (left panel) and 
um/k = 6.4 (right panel). 



VI. TWO OPTICAL MODES 

As mentioned in the introduction, the most studied 
system in optomechanics is the one where an optical cav- 
ity mode is coupled to the motion of a movable end mirror 
PH1J] . The membrane- in-the- middle geometry [15|, [l6| is 
an alternative setup where one avoids having to combine 
a high finesse cavity mirror and a delicate mechanical el- 
ement. In both geometries, it is beneficial to make the 
intracavity photon number and hence the optomechani- 
cal coupling \a\ as large as possible for the observation 
of RPSN. This can be achieved by increasing the power 
of the laser driving the cavity. However, one problem 
with this is that the mean radiation pressure on the me- 
chanical oscillator increases, which leads to an increased 
shift in its equilibrium position. This eventually leads to 
a static instability [lj| [H| where more than one stable 
equilibrium position exist. This effect limits the amount 
of laser power that can be applied. 

It turns out that there is a way to perform the 
membrane-in-the-middle experiment with zero mean ra- 
diation pressure on the membrane. The trick is to drive 
two optical modes in the cavity whose coupling to the 
oscillator is of opposite sign. In other words, the two 
eigenmode frequencies would depend oppositely on the 
oscillator position [l5j]. It should then be possible to in- 
crease the intracavity photon number without worrying 
about the static instability, and thus increasing the possi- 



bility of observing RPSN. Note that this is not possible in 
the movable mirror setup, where the coupling always has 
the same sign. If in addition, both modes are driven by 
the same laser, e.g., by utilizing an acousto-optic mod- 
ulator, the classical noises in the two modes are corre- 
lated whereas the quantum noises are not. The classical 
noises in the two cavity modes will then try to force the 
mechanical oscillator in opposite directions in a synchro- 
nized way, resulting in a small net displacement. This can 
significantly reduce the oscillator fluctuations induced by 
classical laser noise, and thus improve the chances of ob- 
serving the fluctuations induced by photon shot noise. 
We briefly discuss this setup below. Although it is some- 
what specific to the membrane-in-the-middle geometry, 
it may be of interest in other setups as well (22|, HH • 

We now consider two optical modes a and b both cou- 
pled to the mechanical oscillator, 

H = hwyic^c + h (uJc.a + Az) (a) a — (a' a)) (38) 
+ h (wc,6 - Bz) (tfb - (S+S>) + H Ka + H Kb + # 7 , 

where the coupling constants A and B have the same 
sign. If A{a)a) = B{b'b), the mechanical oscillator's 
equilibrium position is the same as in the absence of 
driving the cavity modes. We again assume that the 
cavity is driven from the left, with frequencies WD,a and 
WD,b for the two modes a and b, and define the detun- 
ings A a = WD,a — wc, a and Ab = — ^c.b- We assume 
that the two cavity modes are well separated in frequency 
compared to their linewidths n a and Kb and to the de- 
tunings A a and A 6 , i.e., \uj c , a - wc,b\ > K a,b, \A a ,b\- 

There are now six input modes, a; n ,i(i) and b- ln ^(t), 
with i = L, R, M. We write them as 

ai„, L (i) = e- 1 ^ 1 (a in (t) + L,l) (39) 

k h L(t) = e'^- bt (fc in (t) + 6x) 

a ln ,i(i) = e-^'i^i , $ = R,M 

hn At) = e^'^M , i = R, M , 

where the quantum noise operators have the same 
properties as before. If we assume that the classical 
noises in the two beams are correlated, the coherent 
state amplitudes are ai n (i) = ao + Sx(t) + i5y(t) and 
h n (t) = a,i n (t)bo/ao. We have assumed that ao and 
bo are real. Also, since we are only interested in the 
noises 5x and Sy around the mechanical frequency, we 
can neglect any effects due to different path lengths for 
the two beams. We write the cavity amplitudes as a 

mean and a fluctuating part, a = e~ lUr '- at (d, + d a (t)j 
and b = e _I " D ' bi (b + db(t) J . Linearization of the equa- 
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tions of motion gives 



(y - «A Q ) d a - IOIZ + V K a,RCa,R (40) 

iA b ) <4 + i/3z + ^K fciR £ b 3 



'Kb 

~2 



— (Sx + iSy) + £&. L 
a 



+ y/ K b.M^,b,M 



with a 
are 



- (~ + wm) c - i (a*d a + adl ~ P*db - Pd\j 
Ad and /3 = Bb. The mean cavity amplitudes 



Kax/2 - iA a 
Kb,h/2 - iAi 



«o 



60. 



(41) 



The condition stated above for a zero mean radiation 
pressure on the oscillator is A\d\ 2 — B\b\ 2 . 

The solutions to the equations of motion are given 
in Appendix [B] as well as exact results for the cross- 
correlation (fT5|) in the case of two optical modes. An 
interesting feature is that when, in addition to ^4|a| 2 = 
B|6| 2 , K ai i = Kb,i, A a = Afc are satisfied, the position 
operator z becomes independent of the classical noise Sx 
and Sy, such that S c i,zM vanishes. The properties of 
the quantum and thermal contributions will be the same 
as discussed in Sections Mil and IIVI Although this re- 
quires fine-tuning, one might be able to get close to these 
conditions and thereby reduce the classical contribution 
significantly. 



that in the general case of nonzero detuning, it also dis- 
appears at a specific choice of quadratures for the cross- 
correlation measurement. We have presented figures of 
merit for when the shot noise contribution can be ex- 
pected to dominate over the thermal noise contribution 
in the vicinity of these idealized situations. The choice 
of parameters in the examples we presented is relevant 
to present-day experiments [151 1 1 61 ] , such that the obser- 
vation of radiation pressure shot noise by this method 
should be within reach. 

The relative importance of classical laser noise versus 
quantum shot noise has also been investigated. We find 
that this can be diminished by the use of an asymmet- 
ric cavity and/or a cavity of very high finesse. In addi- 
tion, we proposed a new setup for the membrane-in-the- 
middle geometry involving two optical modes. This setup 
can overcome the problem of static bistability, and can 
potentially also reduce the unwanted classical radiation 
pressure noise significantly. 
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Appendix A: Mathematical details in the case of one 
optical mode 



VII. CONCLUDING REMARKS 

We have presented a detailed theoretical study of a 
proposed experiment designed to observe radiation pres- 
sure shot noise by detecting a correlation between pho- 
ton shot noise in an optical cavity mode and position 
of a mechanical oscillator. The experiment involves the 
measurement of a cross-correlator of two outgoing optical 
field quadratures from the cavity. We have investigated 
the possibility of detecting radiation pressure shot noise 
by this method. 

The cross-correlation measurement has contributions 
from three noise sources: radiation pressure shot noise, 
classical radiation pressure noise and thermal noise in the 
mechanical oscillator. We have determined how the radi- 
ation pressure shot noise contribution differs qualitatively 
and quantitatively from the other contributions to this 
cross-correlator. As pointed out in Ref. [26| . the contri- 
bution from thermal noise disappears when the drive fre- 
quency equals the cavity resonance frequency. We found 



We begin by writing down the solution to the equations 
of motion (JSJ) in Fourier space [35] : 

5M = ^[v^Um'I-^M + xm'M^M) ( A1 ) 

-2cj m (a*xcMCM + axchHC 1 "^]) 
d[u] = -XcH (iaz[w] - CM) ■ 
We have defined the operator 

CM = v^l (Sx[oj] + i5y[w] + (A2) 

+ VkrCrM + V^m£mM , 

the susceptibilites xcM = [ K /2 — i(oj + A)]^ 1 and 
XmM = [7/2 — i(ui — wm)]~ , and the functions 

N[u] = x m 1 HXm 1 *[-^ + 2^mSH, (A3) 
= -i\a\ 2 (xcH - Xcb^D • 
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In the weak coupling limit, — 2iu)MN~ 1 [uj] can be thought 
of as an effective mechanical susceptibility and as 
the optomechanical self-energy. 

It is instructive to write out some of the terms in Equa- 
tion (|26|) . We begin with the (q,z)-term, 



(SX[uj}SYe[uj'}) q , 
x [ <zM£[u/]) 



**A[w'] 



(i R [uj]z[uj'} 



a\- l ^e ie Y,[w] (A4) 
^[-cy]]T^(zM£K]> J 



KRXC M 



The summations are over i — L, R, M. We see that this 
term consists of correlations we wish to detect, namely 
correlations between the quantum vacuum noise of the 
electromagnetic field and position of the mechanical os- 
cillator. We have defined the function 



A[wl 



<*-V&[-w]) , (A5) 



where <f> is the complex phase of the mean cavity ampli- 
tude a. The term (5X[a;]<5Yg[w']) c i !Z contains correlation 
functions of the kind whereas the last term 

in Equation (|26l) is 



(5X[oj]6Y e [uj']} z 



(A6) 



= y/K^\a\- 2 X[uj}A[uj'}(z[uj]z[uj']) 



We next present the terms in Equation (|27|) . The first 
term, which is the term of interest, is 



EM 
N[w] 

AH 
N[-ui] 



(xcH +xhi-u]) 



(A7) 



The contributions from the classical noise in the drive, 
when expressed by the functions 



B±[w] 
D±[u] 

C b ,d[u] 
become 

-Scl,clM 
<Scl,z[w] 



e-^xcHie^XcH (A8) 
e- 16 (1 - KLXoM) ± e w (1 - KLXcH"]) 
B+[u;]5 + [-u;]C x - B_ [cj]B_ [-u]C Y 

B+[u]D + [-u]C x -B-[u]D-[-u]Gy), 



= -^/klkrCb,d[w] 



Afwl 



N[-lj] 



The last term can be written 

£[w]A[-w] 



N[(j]N[-cj] 



(A10) 



M\ 2 + Ixm 1 M I 2 



-4w m k l Cb,bH + 2w M K IxcHI + |Xc[-w 



The first term in S^^u;] is typically the dominant one, 
which represents the fluctuations of the oscillator due to 
thermal noise from the mechanical bath. 

As stated earlier, in the limit 7, \uj — cjm| "C k, cJm> the 
above expressions simplify. The coefficients appearing in 
Equation ([29]) are 

(J^j - A 2 I sw.6 — kAcos( 



(th) 



r(th) 



^( th ) 



/« . 

— SI 
V2 



sin 9 — A cos 



where 



K (th) = -2A 7 ( Tim + ^ ) A^ q > 



(AH) 
(A12) 

(A13) 



The effective phonon number nyi is given by Equation 
(TT21 . with 



and 



^( q ) _ V«l «r I a 1 2 1 Xc [^m] 1 2 1 XC [~^m] 
( K /2) 2 + A 2 



''opt 



«Ixc[~^m]| 2 + k l C b . b [-ujm} 
4 Akwm I XC [wm] 1 2 1 XC [- wm] 1 2 
The coefficients in Equation (I3U1) are 



"" - Ar( q) 7 w M A (A sin 6»+ |cos< 



(A14) 



(A15) 



R 



3A 2 



(I)' 



sm I 



+ A 3 



(3 



A 2 + 4 cos 



r(q) 



2 V V 2 



i 2 +w 2 



M 



(2" 



sin 9 — A cos 6> 



The coefficients in Equation pip are quite complicated 
and therefore not presented here. 



1. Zero detuning 

(A9) 

For I A| <C k, wm and 9^0, the coefficients in Equation 
(1291) become 



(th) -2^/k l k r kA|q;| 2 7 (n M + 5) 



K«/2) 2 +^ M ] 



sin (9 (A16) 
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and I^/R^' = 2ujm/k. For Sq^uj], we get 



Appendix B: Mathematical details in the case of two 
optical modes 



4» = _ vyM; i M11 , 



[(«/2) 



(A17) 



'Ml 



i?^ q) /(i?2 q) 7) = and /{ q) /(4 q) 7) = 1/2- In this case, 
the coefficients in iSd z [w] are relatively simple and given 
by 



R 



(cl) _ 



4^/klkr, k KL|a| 2 Cjs: sinfl 



[(k/2) 



(A18) 



wm/k, J2 /i?2 = — 2u>m/k and 
-1/2. As expected, the oscillator only 



4 cl) /(4 cl) 7) 
/j cl) /(4 cl) 7) = 

sees the amplitude noise and not the phase noise when 
the frequency of the drive is close to the cavity resonance 
frequency. 

It is worth noting that \I% /R^\ = 2u>u/k, whereas 

|4 q) /(i? 2 q) 7)l = V 2 - Unless k > u M , this means that if 
the classical signal is dominant in the real part of S[u>], it 
should also be dominant in the imaginary part (note also 



I2 /i?2 |). Thus, if a sign change 



that \4 th) /R^ h) \ 
is observed in R[u>], but not in I[uj], one might be able to 
conclude that the sign change is due to photon shot noise 
in the cavity and not classical noise in the drive. Another 
possibility is to add classical noise deliberately, in order 
to determine the source of the sign change in R[ui\. 



1. Finite detuning 



The angle 9 C was defined as the quadrature angle 
where the constant R^ vanishes for a given value of 
A/k. For 89 = 9 - 9 C and \89\ < 1, we find 

that R[ th) = K (th ) ((k/2) 2 + A 2 ) 56. The imaginary part 
does not change significantly for 9 around 9 C . Its value 
at 6 = 9 C is 4 th) = if( th WA. We also quote the co- 
efficients appearing in Equations (|30|) at 9 = 9 C . They 
are 



R[ q) = 



(A19) 



R 



(q) 



= A- (q) A 



r(q) 

1 i — 



A 2 

2 



J1 i\l 



The solution to Equations (|40j) are 

^ = ttti\v^(xm*[-^H^} + x M 1 H?? t H) 

-2a; M (a*XaHCaH + ax^[-w]ClH 

-FxbMCM - Pxt[-"K 1 M)] (bi) 
xhH (z/JzM + C&M) ■ 

We have defined 

CaH = V K a,L (Sx[uj] +i6y[u] +(„,lM) (B2) 



d a [ui] 

d b [uj] 



— (5x[uj} + iSy[uj}) + £b,h{u} 
a 



and 



JV2H 
E 2 [w] 

S [w] 



w] + 2wmE 2 [w] 
= E [u;] + E fc [w] 
= -z|a| 2 ( Xa H-x;[-c]) 
= -i\(5\ 2 { Xb [u]-xi[-u]) , 



(B3) 



where the optical mode susceptibilities are Xa/fcM — 

[K a /b/2 — i(u) + A a /(,)] . As stated in Section IVH the 
mean radiation pressure on the oscillator is zero when 
A\d\ 2 = -B|^| 2 - If in addition n a i = Kb,i and A a = it 
is straightforward to check that z becomes independent 
of the classical noise 8x and 8y. This requires fine-tuning, 
but it shows that it might be possible to make the clas- 
sical contribution small enough to be negligible. 

We now return to the general case and imagine mea- 
suring the cross-correlation between the quadrature fluc- 
tuations 



8X(t) = e-^d aM {t) + e**-<2t R (t) (B4) 



and 



5Y e (t) 



(t) 



' l 6,out,L 



(*) (B5) 



in the same 



way 



as before, with d 



j.out.i 



(*) 



e iWD J t (a j -,out,i(«) 



a j, out, i{t))) and e*^ a = a/|a|. SX rep- 
resents the intensity fluctuations in the output of mode 
a on the right hand side of the cavity. SYg is the fluctua- 
tion in an arbitrary quadrature in the output of mode b 
on the left hand side. It is however not essential whether 
mode a or b is detected. Defining S[u] as in Equation 
it still has four contributions as in Equation (|27l) . 
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The quantum contribution is 

X 



(B6) 



H N 2 [uj] 
\a\At\rM 
\P\N 2 [-u] 



(xaM+x:[-w]) 



where we have used the abbreviation Cb^dJu;] = 
S a , + [ix)]I?6, + [— cj]Cx — B a -[uj\Db-[— uj]Cy and similarly 
for CB a ,E[u] etc. Finally, the last contribution is 



with A 6 [w] = - t |/3| 2 (e <A X6H - e- <A x£[-w]), \ = (4>b - 
9) and e*^ 6 = The terms due to classical noise in 

the drive can be expressed by the functions 



B a , ± [u>] = e-^XaM±e^ X ;[-a;] 

B b , ± [u] = ^(e-^x&Mie^Xbha,]) 
ao 



(B7) 



e (1 — k&,lX&M) 

± e «(i-« 6iLX *[- w ])" 
£±M = y^L|a|-B Q;± [w] - v^KftJl^l^iM. 

They become 

Sd,clM = -■ v /Ko,LKo,rCb ,D|,M, 
5 c l,z[u;] = 2uj M y / K a ,R I ^° Cl^gM 



(B8) 



|a| N 2 [w 
1 AbHj 



^H = -^W |a|PaMWa[ _ H (B9) 

x { 7 ( nth + £) (^mI-MP + Ixm'HI 2 ) 

H-^Cb.bM +2^[k |q;| 2 (|XaM| 2 + |Xa[-w]| 2 ) 

+« fc |/3| 2 (|x h H| 2 + |x 6 h^]| 2 )' 



Again, we point out that when A\a\ 2 — B\b\ 2 , A a — 
/t 0i j = /tfe,,, the functions -E±[w] are identically zero and 
hence 5 c i, z M = 0. In this case, when analyzing the other 
contributions, one finds the same results as in the case 
with one optical mode, only with modified 7 op t and Sojm- 
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